Quantum vs. classical hyperfine-induced dynamics in a quantum dot 
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In this article we analyze spin dynamics for electrons confined to semiconductor quantum dots 
due to the contact hyperfine interaction. We compare mean-field (classical) evolution of an electron 
spin in the presence of a nuclear field with the exact quantum evolution for the special case of 
uniform hyperfine coupling constants. We find that (in this special case) the zero-magnetic-field 
dynamics due to the mean-field approximation and quantum evolution are similar. However, in a 
finite magnetic field, the quantum and classical solutions agree only up to a certain time scale t < r c , 
after which they differ markedly. 
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I. INTRODUCTION 

Prospects for future quantum information processing 
with quantum-dot-confined electron spins* have encour- 
aged a series of recent experimental efforts. These ef- 
forts have resulted in several very significant achieve- 
ments, including single-electron confinement in vertical 2 - 
and lateral single 3 - and double 4 ^ gated quantum dots, 
the demonstration of spin-dependent transport in dou- 
ble dots^^ and exciting effects arising from the contact 
hyperfine interaction with nuclear spins in the host ma- 
terial, including coherent undriven oscillations in spin- 
dependent transport^ lifting of the spin-blockade^ en- 
hancement of the nuclear spin decay rate near sequential- 
tunneling peaks)2ii£ and notably, decay of coherent os- 
cillations between singlet and triplet states as well as 
the demonstration of two-qubit gates in double quantum 
dots . 11 ' 12 Very recently, the hyperfine interaction has also 
been identified as the source of decay for driven single- 
spin Rabi oscillations in quantum dots . 13 ' 14 

In spite of rapid progress, there are still many obstacles 
to quantum computing with quantum dots. In particular, 
the inevitable loss of qubit coherence due to fluctuations 
in the environment is acceptable in a quantum computer 
only if the error rates due to this loss are kept below 
10~ 3 — 1CP 4 errors per operation.— This requirement is 
particularly difficult to achieve since it means that inter- 
actions must be strong while switching so that operations 
can be performed rapidly, but still very weak in the idle 
state, to preserve coherence. 

For an electron spin confined to 
decoherence can proceed through 
the electromagnetic environment 
interactio n 16 ' 17 ' 18 ' 19 ' 20 or through 
interaction with nuclei in 

material, which has been siiown extensively m 
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and experiment . 1 1 ' 12 i 44 ] 45 i 46 i 47 i 48 Due to the primarily 
p-type nature of the valence band in GaAs, hole spins 
(unlike electron spins) do not couple to the nuclear 
spin environment via the contact hyperfine interaction, 
although they can still undergo decay due to spin-orbit 
coupling. The decay may still occur on an even longer 
time scale than for electrons^ 9 - which suggests the 
dot-confined hole spin may be another good candi- 
date for quantum computing. Alternatively, quantum 
dots fabricated in isotopically purified 28 Si^ or 12 C 
nanotube o 51 ' 52 ' 53 would be free of nuclei with spin, and 
therefore free of hyperfine-induced decoherence. 

While the field of quantum-dot spin decoherence has 
been very active in the last few years, there still re- 
main significant misconceptions regarding the nature of 
the most relevant (hyperfine) coupling, particularly, the 
range of validity of semiclassical spin models and tradi- 
tional decoherence methods involving ensemble averaging 
have been called into question for a single isolated quan- 
tum dot with a potentially controllable environment. We 
address these issues in section [TT1 



II. HYPERFINE INTERACTION: QUANTUM 
AND CLASSICAL DYNAMICS 



Exponential decay of the longitudinal and transverse 
components of spin is typically measured by the decay 
time scales T\ and T2, respectively^ The longitudinal 
spin relaxation rate 1/Ti due to spin-orbit interaction 
and phonon emission is significantly reduced in quantum 
dots relative to the bulk in the presence of a weak Zee- 
man splitting B = |B| and large orbital level spacing 
hui (l/Ti oc B 5 /(hui a ) 4 )^.^ This decay time has been 
shown to be on the order of T\ ~ 1 ms in gated GaAs 
quantum dots at B sa 8 T^ and to reach a value as large 
^as Ti = 170 ms at low magnetic fields (B = 1.75 T)^ 6 . 
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Furthermore, since dephasing is absent at leading order 
for fluctuations that couple through the spin-orbit in- 
teraction, the T2 time due to this mechanism is limited 
by the T\ time (T2 = 2Ti)i& (we note that corrections at 
higher order in the spin-orbit interaction can lead to pure 
dephasing, although these corrections are only relevant at 
very low magnetic field a 57 ' 58 ). Unlike the spin-orbit in- 
teraction, the hypcrfinc interaction can lead to pure de- 
phasing of electron spin states at leading order, resulting 
in a relatively very short decoherence time r c « 1 — 10 ns 
due to non-exponential (Gaussian) decay 23 ' 26 To perform 
quantum-dot computations, this and any additional de- 
cay must be fully understood and reduced, if possible. 

The Hamiltonian for an electron spin S in the lowest 
orbital level of a quantum dot containing nuclear spins is 

# hf = S-(B + h); h = J]A l I l , (1) 

i 

where A, = Avq \ipo(ri)\ 2 is the contact hyperfme cou- 
pling constant to the nuclear spin at site i, vq is the 
volume of a crystal unit cell containing one nuclear 
spin, and A w 90 fj,eV is the weighted average hyper- 
fine coupling constant in GaAs, averaged over the cou- 
pling constants for the three naturally occurring radioiso- 
topes 69 Ga, 71 Ga, and 75 As (weighted by their natural 
abundances),— all with total nuclear spin / = 3/2. The 
nuclear field in H^f is given by the quantum "Over- 
hauser operator" h. Although an exact Bethc Ansatz 
solution exists for TJhf,— using this solution to per- 
form calculations for the full coupled quantum system 
of TV « 10 4 — 10 6 nuclei and one electron in a quan- 
tum dot can be prohibitively difficult Since the Over- 
hauser operator h is a sum of a large number N of spin-/ 
operators, one expects that under certain conditions its 
quantum fluctuations can be neglected and the opera- 
tor h can be replaced with a classical Overhauser field 

h — > E iv S - 11 - 23 - 26 - 33 - 37 - 38 - 44 - 45 - 61 - 62 ^. 64. 65. 66. 67. 68 JJ 0W _ 

ever, this approximation can accurately describe the 
electron-spin dynamics only at times t < r c , where 
t c = N v I A and 77 > O 6 ^^ after which effects of quantum 
fluctuations of the Overhauser operator set in. The nu- 
clei in GaAs are indeed quantum objects, which could be 
verified, in principle, by demonstrating that they can be 
entangled, as is done in spin-state squeezing experiments 
that have been performed on atomic ensembles.— The re- 
placement h — > B^r is therefore not exact and there are 
several cases in which the electron-spin dynamics at times 
t > t c differ markedly for quantum and classical nuclear 
fields. In particular, without performing an ensemble 
average over initial Overhauser fields, the classical- field 
picture predicts no decay of the electron spin. This is 
in direct contradiction to analytica l 31 ' 41 ! 42 ' 70 and exact 
numerica l ' 35 studies that show the quantum nature of 
the nuclei can lead to complete decay of the transverse 
electron spin, even in the presence of a static environ- 
ment (fixed initial conditions). Additionally, quantum 
"flip-flop" processes can lead to dynamics and decay of 
the electron spin in the quantum problem, even for initial 



conditions (e.g., a fully-polarized nuclear spin system) 
that correspond to a fixed-point of the classical equa- 
tions of motion . 23 ' 31 ' 36 In fact, it can be shown that any 
decay of the electron spin for pure-state inital conditions 
will result in quantum entanglement between the electron 
and nuclear spin systems . 27 ' 28 This entanglement has re- 
cently been highlighted as a source of spin-echo enve- 
lope decay in the presence of the hyperfme interaction^ 3 - 
Finally, even the ensemble-averaged standard classical 
(mean-field) electron-spin dynamics show large quanti- 
tative differences relative to the exact quantum dynam- 
ics at times t > t c and in a very weak magnetic field, 
although an alternative mean-field theory involving the 
P-representation for the density matrix shows promise.— 

While the classical and quantum dynamics diverge in 
many cases, the classical-field replacement h — > B^v will 
be valid up to some time scale, providing a range of 
validity for the classical dynamics. In this article, we 
aim to shed light on this range of validity of the classi- 
cal solution. As a test of the classical-dynamics picture, 
we can compare quantum and classical dynamics of an 
electron spin in the simple case of uniform coupling con- 
stants Ai = 7. When the coupling constants arc uniform, 
an exact solution to the quantum dynamics (see Rcfs. 
241(721 for the |B| = case) can be evaluated and used 
to compare with an integration of the equivalent classi- 
cal equations of motion. For uniform coupling constants, 
the nuclear Overhauser operator from Eq. JT]) becomes 
h = 7K, where Ai = 7 = A/N and K = £\ is tlic 
collective total spin operator for N 3> 1 nuclear spins. 

The initial state of the system is taken to be an arbi- 
trary product state of the electron and nuclear system: 

|V(0)> = |Vs(0)) ® |^jr(0)> , (2) 

K 

= Y, (al\V,K,m) + al\l;K,m)), (3) 

m=-if 

where \a]K, m) is a simultaneous eigenstate of S z , K z 
(we take the direction of the external field B to de- 
fine the z-axis), and K • K (with eigenvalues ±1/2 for 
(7 =|,|, to, and K(K + 1), respectively). For compari- 
son with the classical spin dynamics, we choose the col- 
lective nuclear spin to be initially described by a spin 
coherent state, given by |V*(0)) = e~ iK ^ K \K, K) = 

Em^St* 9 *") \K,m), where d { £ K (e K ) is the Wigner ro- 
tation matrix^ 3 - and the electron spin is in an arbitrary 
initial state |Vs(0)> = cos((9 s /2) ||) + e^ s sin(0 s /2) ||). 
The initial conditions are then completely determined by 
the three angles 0s,<f>s, and 9k- These initial conditions 
allow for an arbitrary relative orientation of the spin and 
magnetic-field vectors, since the azimuthal angle for K 
(4>k) can be set to zero with an appropriate shift in (jjs- 
(f)' K = 0, 4>' s = 4>s ~ 4>k- At any later time t, the wave 
function is given by 

K 

\ij(t))= E ( a m(<) It; K, m) + ai(t) ||; K, m» . 

rn=-K 

(4) 
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From the time-dependent Schrodingcr equation 
id t \ip(t)) = Hhf\ip(t)) (setting h = 1), we find the 
set of coupled differential equations determining the 
coefficients {a m (t), aj^(t)} . For m = — K...K — 1, 



„T _ 



--(B + -/m)al 



■1 n - i 

1 2 U Km+l a m+V 



Kn+i = ^{B + lim + l^a 



I 

m+l 



I 

V 



(5) 
, (6) 



The solutions to Eqs. ([5]), ([6]) and the expressions in 
Eqs. ©, © for the coefficients {a4(*), a^(t) : m = 
~ K . . . K} constitute a complete exact solution for the 
dynamics of the wave function \tf>(t)} at any later time t > 
0. We solve Eqs. ([5]) and ([6]) by Laplace transformation 
to obtain 



where C 



Km 



(Km±l\K± \Km) 



y/K{K + 1) — m(m ± 1). These equations arc sup- 
plemented by two equations for the stationary states 
|t; K, K) and \i;K, —K) with dynamics: 



a ] K {t) = expi-J(S + 7 ^)iUt(0) 



(7) 



at K (t) = e X p\-(B- 7 K)tjat K (0). (8) 
I 



,+i(*) 



LdKr 



e l 2* |a4(0) cas{u Km t) - i fa^(0) 
e*** \ a l m+1 (0) cos(aW) + i ( a l m+1 (0) 



B + 7 ( m+ _) 



+ «r i „+l(°)7^ m+1 



s + 7 ( m + _) 



(B + 7m) (B + 7(01 + 1)) + 7 C 



y Km-\-l Km 



-ai(0) 7 C+ f 

1/2 



sin(cj g mQ 
sin(cJK m f) 



(9) 
(10) 

(11) 



With the coefficients {a£,(t), a^i) 



-if... a:} 



in hand, we can evaluate the expectation values of all 
spin components exactly: (S) t = (tp(t)\ S \ip(t)) , (K) t = 

(m\K\m)- 

To evaluate the classical spin dynamics, we perform a 
mean-field decomposition of the Hamiltonian given in Eq. 
(fTJ by rewriting the spin operators as S = (S) t + SS and 
K = (K) f + (5K. We then neglect the term that is bilinear 
in the spin fluctuations (oc SS ■ (5K) and approximate the 
spin expectation values by their self-consistent mean-field 
dynamics (S) t ~ s(<), (K) t w k(t), where s and k are 
classical time-dependent vectors of fixed length^ Up to 
a c- number shift, this gives the (time-dependent) mean- 
field Hamiltonian 

H ml (t) = (B + 7 k(t)) • S + 7 s(i) • K. (12) 

The mean-field dynamics are now given by the Heisen- 
berg equations of motion for the spin operators: S = 
i[HjDf(t),S], K = i [H m f (t), K], with the replacements 
(S> t «s(i), (K> t «k(i): 

s(t) = (B+ 7 k(i)) xs(t), (13) 

k(i) = - 7 k(t) x s(i). (14) 

An exact analytical solution to Eqs. (|13I14[) is known^ 
However, instead of repeating this solution here, we solve 



Eqs. (|13I14[) by numerical integration for direct compar- 
ison with the exact results given above. The mean-field 
and quantum dynamics are shown in Fig. [Tj for four val- 
ues of the Zeeman splitting B = |B|. We compare the 
two solutions using the correlation function 



C(t) 



1 
T 



t+T 



dt'- 



2{S x ) t , s x (t') 



(S x ) 



A2 



(15) 



where we average over the time interval T = O.l/i/7 to 
remove rapid oscillations. C(t) = 1 if the exact solu- 
tion and mean-field approximation are identical (s x (t) = 
(S x ) t ) over the time interval (t,t + T). C(t) < 1 indicates 
that the two solutions differ. While the zero-magnetic- 
field dynamics appear to be well reproduced by the mean- 
field approximation, at least at short time scales, the 
high-field solution decays rapidly, which can not appear 
in the classical dynamics unless averaging is performed 
over the initial conditions^ There is a partial recurrence 
of the correlator at a time scale given by the inverse level 
spacing for the quantum problem r p = 2irh/^/, but the 
recurrence is only partial since at this time the quantum 
and classical solutions have already gone out of phase. 

It is relatively straightforward to understand the differ- 
ence in the high-field and low-field behavior shown in Fig. 
[TJ At zero magnetic field, the total spin J • J ( J = K + S) 
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FIG. 1: (Color online) (a) Correlation between the mean-field and exact quantum solution C(i) (defined in Eq. (|15[) : C(t) = 1 
indicates perfect agreement between the mean-field and quantum solutions) for evolution of an electron spin in the presence 
of a total bath spin K = 50 and magnetic field B = (black dash-dotted line, showing the weakest decay), B — 57 (blue 
dotted line), B = IO7 (green dashed line) and large magnetic field B — 5OO7 (red solid line, showing rapid decay). The inital 
conditions were 0$ = 7r/2, <j>s = 0, 6jc = 0.3-7T (see the discussion following Eq. @). We also show the exact quantum evolution 
{S x ) t (solid line) and mean-field approximation s x {t) (dashed line) for (b) B = and (c) B — IO7. 



commutes with the Hamiltonian, so if the nuclear spin 
system begins in an eigenstate of K • K, only a single 
frequency exists in the quantum dynamics, correspond- 
ing to the difference in energies with J = K ± l/2i 24 i 28 
Thus, in this case the quantum dynamics corresponds 
to simple periodic precession, and mimics the classical 
dynamics for K >> 1 (sec Fig. QJb)). However, the 
states of fixed J are manifold degenerate. If a term is 
added to the Hamiltonian which does not commute with 
J • J (in this case, the electron Zeeman term BS Z ), many 
more frequencies are involved in the quantum dynamics, 
which can lead to decay in the quantum solution, while 
the classical solution continues to describe simple elec- 
tron spin precession (see Fig. [ljc)). In a large magnetic 
field (B ^> 7|Kj_|), it is straightforward to evaluate the 
decay in the quantum mechanical solutio n 31 ) 76 



(S+) t « (S + ) fflcp^-^+i [5 + 008(^)7^*^6) 



decay time 



7 V K [I -cos 2 (6 K )] 



(17) 



N 



A 



(18) 



Since the classical dynamics at times t < t c describe 
simple precession for fixed initial conditions, any decay 
in the quantum solution signifies a disagreement between 
the quantum and classical problems. Thus, the mean- 
field solution will give an accurate description of the full 
quantum dynamics only for times t < r c , with r c given 
by Eq. HI]). 

The crossover from precession to decay of the quan- 
tum solution with the addition of a magnetic field sug- 
gests that the uniform coupling-constants picture should 
only be used with caution, since the Hamiltonian in Eq. 
([I]) also does not commute with J • J when the coupling 
constants vary from one nuclear-spin site to the next 
(as is true in a quantum dot). Indeed, in the presence 
of randomly-varying coupling constants, the straightfor- 
ward mean-field electron-spin dynamics at times t > r c 
are quantitatively very different from the exact quantum 
dynamics at weak magnetic fields B — > Oi^I 



The x-component of spin is then given by the real part 
(S x ) t = Rc[(5+)J. We consider the hyperfine problem 
with / = 1/2. When the initial nuclear-spin coherent 
state is generated by rotating the spins from a fully- 
polarized state such that K is maximal (as in Ref. [3lh . 
we then have K = N/2. In addition, 7 = A/N and for 
nuclear spin polarization p — cos(#x) *C 1 this gives the 



III. CONCLUSIONS 

We have presented an exact solution for the prob- 
lem of an electron spin interacting with a large bath 
of spins with uniform Heisenberg coupling. This exact 
solution has been compared to the corresponding mean- 
field (classical spin) model. We have seen that the mean- 
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field and quantum solutions show striking agreement at 
times shorter than the transverse-spin correlation time 
r c , which diverges at zero magnetic field. This diver- 
gence, however, may only be due to the assumption of 
uniform coupling constants, which is unphysical for a 
quantum dot with strong confinement. 

In this work we have focused on a comparison of dy- 
namics for fixed initial conditions of the quantum and 
classical problem. Some of the quantum behavior, in- 
cluding Gaussian decay, can be recovered with an aver- 
age over classical solutions . 23 ' 26 An intriguing question 
therefore remains: How much of the quantum dynamics 
can be obtained by averaging over classical solutions with 
different initial conditions? 



Note added: Recently, a related preprinUi has ap- 
peared in which the authors use the exact solution for 
uniform coupling constants to evaluate the z-component 
of electron spin, complementing earlier predictions for the 
bath-polarization dependence of decoherence in single 3 ^ 
and double dots4^ 
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(Color online) (a) Correlation between the 
mean-field and exact quantum solution 
C{t) (defined in Eq. (JTSJ); C{t) = 1 
indicates perfect agreement between the 
mean-field and quantum solutions) for 
evolution of an electron spin in the pres- 
ence of a total bath spin K = 50 and mag- 
netic field B = (black dash-dotted line, 
showing the weakest decay), B = 57 (blue 
dotted line), B = IO7 (green dashed line) 
and large magnetic field B = 5OO7 (red 
solid line, showing rapid decay). The ini- 
tal conditions were 9s = 7r/2, <ps = 0, 
9 K = 0.37T (see the discussion following 
Eq. ([3])). We also show the exact quantum 
evolution (S x ) t (solid line) and mean-field 
approximation s x (t) (dashed line) for (b) 
B = and (c) B = IO7 



